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Abstract. The light-by-light contribution to the anomalous magnetic moment of muon (g — 2) M from 
the hadronic exchanges in the neutral pseudoscalar meson channel is calculated in the nonlocal chiral 
quark model. The full kinematic dependence of the meson-two-photon vertices from the virtualities of 
the mesons and photons is taken into account. The status of various phenomenological and QCD short- 
distance constraints is discussed and the comparison with the predictions of other models is performed. It 
is demonstrated that the effect of the full kinematic dependence in the meson-photon vertices is to reduce 



X 



the contribution of pseudoscalar exchages to ajj 
most of previous estimates. 



1 Introduction 

The study of the anomalous magnetic moments (AMM) of 
leptons, a — (g — 2)/2, has played an important role in the 
development of the standard model. At the present level 
of accuracy, the muon AMM, a M , leads to the sensitivity 
to effects of new physics, with particles in the mass range 
of a few hundred GeV/c 2 . 

Comparison of the experimental measurements of the 
muon AMM pQ with the predictions of the standard model 
(see, for example [2]) shows deviation by approximately 
3ct. The uncertainties of the standard model value for 
are dominated by the uncertainties of the hadronic contri- 
butions a?, trong , since their evaluation involves methods of 
nonpcrturbative QCD at large distances. The contribution 
to the muon AMM to leading order in the QED coupling 
constant a comes from the hadronic vacuum polarization 
and the next-to- leading corrections consist of the contribu- 
tions which are the iteration of the leading order term plus 
the contribution from the light-by-light (LbL) scattering 
process. The absolute values of the leading and next-to- 
leading terms are differed by one order of magnitude, but 
the theoretical accuracy of their extraction is comparable 
and dominates the overall theoretical error of the standard 
model result. The hadronic vacuum polarization contri- 
bution to the muon AMM is known with accuracy better 
than one percent owing to the phenomenological analysis 
of the inclusive e + e~ — > hadrons and r — > hadrons data 
[5]. This, however, is not the case for the LbL contribution 
known from various QCD motivated approaches with ac- 
curacy of order 50%. It is the LbL contribution, which will 
dominate the theoretical error of the standard model pre- 



by approximately factor 1.5 in comparison with the 



diction in the near future, and more realistic approaches 
are needed for a better understanding. 

In general, the LbL scattering amplitude is a compli- 
cated object for calculations. It is a sum of different dia- 
grams, the quark loop, the meson exchanges, the meson 
loops and the iterations of these processes. Fortunately, 
already in the first papers [H03E] devoted to the calcula- 
tion of the hadronic LbL contribution, it was recognized 
that these numerous terms show a hierarchy. This is re- 
lated to existence of two small parameters: the inverse 
number of colors 1/N C and the ratio of the characteris- 
tic internal momentum to the chiral symmetry parame- 
ter m jU /(47T /„■) ~ 0.1. The latter suppresses the multiloop 
contributions, so that the leading contributions are due to 
the quark loop diagram and the two-loop diagrams with 
mesons in the intermediate state. Here the contribution of 
the diagram with intermediate pion is enhanced by small 
pion mass in the meson propagator. 

Different approaches to the calculation of the pseu- 
doscalar meson contributions to the muon AMM from 
the light-by-light scattering process are used. These ap- 
proaches can be separated in two groups. The first one 
consists of various versions of the vector meson domi- 
nance model (VMD) [tJllTllgMTU] and the second one of 
QCD inspired effective models (EM). They include differ- 
ent versions of the Nambu- Jona-Lasinio model [TlTT!?] , the 
models based on nonpcrturbative quark-gluon dynamics, 
like the instanton liquid model [T3], the Schwinger-Dyson 
model [HUH (DSE) or the holographic models (HM) [15]. 
To reduce the model dependence of various approaches, 
different constraints on their parameter space are em- 
ployed. One kind of important constraints on the models 
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follows from the phenomenology of the two-photon widths 
of the pseudoscalar mesons r (P — » 77) and their transi- 
tion form factors Fp 77 . (— Mj,;0,q 2 ) first emphasized in 
[6 . Another set of constraints follows from the large mo- 
mentum asymptotics for the meson transition form factors 
07] and for the total light-by-light scattering amplitude 
considered in [81H3]. obtained using perturbative QCD. In 
addition, the model amplitudes have to be consistent with 
the 4-momentum conservation law. Finally, the model cal- 
culations should be tested by calculations of the hadronic 
vacuum polarization to the muon AMM |17U18U19U20[|21j . 

The present work is devoted to the calculation of the 
pseudoscalar mesons contribution of the hadronic light- 
by-light scattering process to the muon AMM within the 
nonlocal chiral quark model (NxQM) [22 , 23 ,24 . In ear- 
lier work [13) the contribution of the pion exchange was 
estimated within this approach in the chiral limit. Now we 
include the effect of the current quark masses and extend 
the calculations also to 77 and 77' meson exchanges. As was 
emphasized in |13j . one of the main advantage of the ef- 
fective model approaches like NxQM is to control the full 
kinematical dependence of the LbL amplitude. The paper 
is organized as follows. In Sect. 2 the basic elements of the 
effective model is given. In Sect. 3 the properties of the 
three-point amplitude with arbitrary pseudoscalar meson 
and photon virtualities are considered. The numerical re- 
sults on the light-by-light contribution in the pseudoscalar 
channel obtained within the nonlocal chiral quark model 
are given in Sect. 4. Sect. 5 is devoted to discussion of the 
constraints for the LbL amplitude that allow to dimin- 
ish the model dependence of the estimates. Sect. 6 and 7 
contain our results in comparison with other model calcu- 
lations and conclusions. 



where A a are the Gell-Mann matrices for a = 1,..,8 and 
Ao = \f2jil. The nonlocal structure of the model is in- 
troduced via the nonlocal quark currents 

J m( x ) = / d 4 x 1 d i x 2 f(xi)f(x 2 ) q{x - x{) r^q(x + x 2 ), 



where M = S,P and — A°, Jp = ij 5 \ a , and f(x) is a 
form factor reflecting the nonlocal properties of the QCD 
vacuum. 

The model fTJ can be bosonized using the station- 
ary phase approximation which leads to the system of 
gap equations for the dynamical quark masses m^j (i — 
u, d, s) 



m d , u + GS U + —S U S B = 0, 
m diS + GS S + —Si = 0, 



with 



-8N C 



d%k 



f 2 (k 2 )m t (k 2 ) 



{2ttY A(fc 2 ) 



(2) 



(3) 



where m^fc 2 ) = m Cti + md } if 2 (k 2 ), Di{k 2 ) = k 2 + m 2 (k 2 ) 
is the dynamical quark propagator obtained by solving 
the Schwinger-Dyson equation, f(k 2 ) is the nonlocal form 
factor in the momentum representation. 

The vertex functions and the meson masses can be 
found from the Bethe-Salpeter equation. For the separable 
interaction (JTJ) the quark- antiquark scattering matrix in 
pseudoscalar channel becomes 



T = t(p 2 )<5 4 (p! + P2 - (ps+pt)) n f(pf), 



2 N^QM Lagrangian, T matrix and r\ — rj 
mixing 

The Lagrangian of the SU(3) x SU (3) chiral quark model 
has the fornQ 

C = q(x)(id - m c )q(x) + ^[J a s (x)J a s {x) + J a P {x)J a P {x)\ 

- jT abc [r s (x)J h s (x)J c s (x) ~ iJ a s {x)J P {x)J P {x)], (1) 

where q (x) are the quark fields, m c is the diagonal ma- 
trix of the quark current masseaj, G and H are the four- 
and six-quark coupling constants. Second line in the La- 
grangian represents the Kobayashi-Maskawa-t'Hooft de- 
terminant vertex with the structural constant 



1 For simplicity in this work we do not consider an extended 
model that includes other structures besides the pseudoscalar 
(P) and scalar (S) ones. 

2 We consider the isospin limit m CiU — m c> d 7^ m Cl a- 



T(P 2 ) = «7sA fc 



-G 



n( P 2 ; 



(4) 



where pi are the momenta of external quark lines, G and 
II (p 2 ) are the corresponding matrices of the four-quark 
coupling constants and the polarization operators of pseu- 
doscalar mesons (p = pi+p 2 = P3+Pi)- The meson masses 
can be found from the zeros of determinant det(G _1 — 
II(— M 2 )) = and the T-matrix for the system of pseu- 
doscalar mesons can be expressed in the form 



T(P 2 )= E 



Va(p 2 )®Vg(p 2 ) 

-(p 2 + M 2 ) ' 



(5) 



where M a are the meson masses, V a (p 2 ) are the vertex 
functions (V a (p 2 ) = 7°V^(p 2 )7°). In general case of three 
unequal quark masses it is necessary to solve the tt° — 77 — 77' 
(or 77 — 77 s — 77 s ) system factorizing in the isospin limit 
into the ir° and 77 — 77' systems. For the 77 — 77' system 
it is convenient to diagonalize the scattering matrix by 
orthogonal transformation 



77 \ _ / cos 9 — sin 9\ ( rjs 
rj' J \sin6> cos£> J 1 770 



(6) 
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As a result, the mesonic vertex functions introduced in ([5]) 
become 

V> ( P 2 ) = i^{p 2 )X 3 , (7) 
V v (p 2 ) = il59 v (p 2 ) (A 8 cos6»(p 2 ) - A sm9(p 2 )) , 
V v > (p 2 ) = iTodv'iP 2 ) ( A s sm9{p 2 ) + A cos6{p 2 )) , 

where g a (p 2 ) and #(p 2 ) are the meson renormalization con- 
stants and mixing angles depending on the meson virtu- 
ally. The renormalization constants arc defined through 
the unrenormalized meson propagators D a (p 2 ) as 



Therefore 9 r , 



9 2 a (p 2 ) = -(p 2 + M 2 a )D a (p 2 ). 
The G-matrix elements have the form 



(8) 







H 


Goo 


= G- 


J 






H 


Gs8 


= G + 


¥ 


Go8 


= Gso 






= G33 


= ( 



(9) 



V2 



H(S U — S s ), 



H 



The matrix II (p 2 ) is diagonal in the quark- flavor basis and 
in the singlct-triplet-octet basis it is given by 



n 00 ( P 2 ) = 1 {2n uu (p 2 ) + n ss { P 2 )) , 
n 88 ( P 2 ) = 1 (n uu (p 2 ) + 2n ss ( P 2 )) , 

n 08 ( P 2 ) = n 8a ( P 2 ) = ^ (n uu ( P 2 ) - n ss ( P 2 )) 
n,{p 2 ) = n 33 ( P 2 ) = n uu (p 2 ), 

where the polarization operators are 

d%k f 2 (k 2 + )f 2 (k 2 _) 



(10) 



(2tt)4 Di(k 2 + )Dj(k 2 _)- 
x [(k+ ■ k-) + m i (k\)m j (k 2 _)] 



(11) 



and k± — k±p/2. The unrenormalized mesonic propaga- 
tors for the pseudoscalar mesons are 

D^(p 2 ) = -G^+II 7r {p 2 ), 
1 



A = -G 88 /det(G) + n QQ (p 2 ), 
B = +G 08 /det(G) + n 0S (p 2 ), 
C = -G o/det(G)+n 88 (p 2 ), 
dct(G) = GooGss — Gq 8 . 

The meson mixing angle depends on the virtuality 



(12) 



8(p 2 ) — — arctan 



2B 



A-C 



;0(A-C). 



(13) 



?(-M 2 ) and e v > = 0(-M 2 ) are different 
for the on-shell r\ and rj' mesons. 

For numerical estimations we use the Gaussian nonlo- 
cal form factor 



f(k 2 ) = eM~k 2 /2A 2 ), 



(14) 



and the model parameters obtained in |24j . In that work 
these parameters (the current quark masses m c ,i, the cou- 
pling constants G and H, and the nonlocality scale A) are 
fixed by requiring that the model reproduces correctly the 
measured values of the pion and kaon masses, the pion de- 
cay constant / w , and the rj' mass (parameter sets G/, Gjv) 
or the T]' — > 77 decay constant gy 77 (sets G/j, Gju). For 
completeness, we present the parameter sets Gi-iv and 
the basic meson properties in Tables Q] and [5] as they were 
determined in [53]. The sets G/, G/y vary by different 
input for the nonstrange current quark mass, while G77, 
Gin are two solutions of the same fitting procedure. 



3 Anomalous triangles 

By using the nonlocal chiral quark modejf], the triangular 
diagram with external pseudoscalar meson and two pho- 
ton legs with arbitrary virtualities (Fig. [1} can be written 

as 

A(7*(9i,ex) 7 *(<?2,e 2 )^P*(p))= (15) 

■ 2 U V P (7 T7l ( 2 2 2\ 

-te e ll „ pa e^e 2 q'[q 2 Fp^, T (p ;q 1 ,q 2 ) , 

with the photon momenta q± t2 and the polarization vectors 
6^2, p = qi + <?2 - For different pseudoscalar meson states 
one has 

„2. 2 „2\ 



F Kn*r (p 2 '1i^ ( 1 2 ) = 9tt(p 2 )F u (p 2 ;ql,q 2 ) 

I 2 2 2\ 9v(P 2 ) 



3V3 

(5F U { P 2 ;qlq 2 2 ) ~ 2F S ( P 2 ;qlq 2 2 )) cos9(p 2 )- 

- V2 (5F U (p 2 ; q 2 , q 2 ) + F s (p 2 ; q 2 ,q 2 )) sm9(p 2 ) 

F - , (r?-a 2 a 2 ) - i^E^l x 
r^/« 7 » 7 . \ja ,qx,q 2 ) — x 

x (5F U { P 2 ;qlq 2 2 ) - 2F S ( P 2 ;qlq 2 2 )) sm9(p 2 )+ 

+ V2 (5F U (p 2 ; ql q 2 2 ) + F s (p 2 ; q\,ql)) cos8(p 2 ) 



(16) 



with 



3 For the pion, the corresponding vertex in the chiral limit 
has been considered in [25] . 
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set 




md,u 


m c , s 


m d , s 


A 


GA 2 


HA'° 




[MeV] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 






G X 


8.5 


304.5 


223 


427 


709 


21.986 


-1670.19 


Gii 


8.5 


304.5 


223 


439 


709 


22.898 


-1557.28 


Gin 


8.5 


304.5 


223 


422 


709 


21.605 


-1717.59 


Giv 


7.5 


287.5 


199 


408 


768 


20.896 


-1721.69 



Table 1. The model parameter sets obtained in [24| . 



set 


M w 
[MeV] 


M v 
[MeV] 


M„, 
[MeV] 


5^77 

[GeV" 1 ] 


9vn 
[GeV" 1 ] 


9r) '77 

[GeV" 1 ] 


Gi 


138.9 


516.5 


958.4 


0.2706 


0.3082 


0.3752 


Gii 


138.9 


505.4 


878.6 


0.2706 


0.3259 


0.3401 


Gin 


138.9 


520.7 


1006.4 


0.2706 


0.3011 


0.3489 


Giv 


139.0 


522.1 


> 739.7 


0.2713 


0.3068 




exp 


134.9766 


547.8533 


957.78 


0.2744 


0.2726 


0.3423 




±0.0006 


±0.024 


±0.06 


+0.009 
-0.008 


±0.008 


±0.014 



Table 2. The basic meson properties for different parametrizations obtained in |24] , 



/(*?)/(*§) 



Ji 



A(*i)A(A3)A(* 2 

m i (k 2 )-mf ) {k x ,k)J x - mf\k 2 ,k)J 2 
.„2 M gl(fcgi)(fcigi) - qj(kq2)(k 1 q 2 ) 

„2„2 



QlQ2 - (9l<?2) 2 



(17) 



.„2 j_ ql{kq 2 ){k 2 q 2 ) - ql{kqi){k 2 qi) 
„2„2 



QlQ2 ~ (lll2) 



where k\ = k ± q x , k 2 = k — q 2 , m^(k,p) — (rrii(k 2 ) — 
rrii(p 2 j)/(k 2 —p 2 ) is the first order finite-difference deriva- 
tive. From (fTT|) one can easily obtain the expressions for 
some special kinematics 



F^qlqlO) 



d%k f(k 2 )f(k 2 ) . 
(2tt)4 Dl {k 2 )D 2 {k 2 ) 



(18) 



m,(l 2 ) - mf '[ku k)Ji - m' t (k 2 )J 2 



Ji (k,qi 

^(0;0,0) = 
— 8m r 



k 2 - 
1 



(kqi)+ 3 
(kqif 



<7? 
1 

47T 2 



(19) 



d%k mi(k 2 ) - 2m' i (k 2 )k 2 



(27T) 



Dim 



In particular, the kinematics displayed in (1181) is of special 
interest for the hadronic exchange LbL calculations (see 
Eq. (|!20|) below). In (|19p the first term is due to the axial 
anomaly, while the second term represents the correction 
due to explicit breaking of the chiral symmetry by current 
quark mass. 



4 Light-by-light hadronic contribution to the 
muon AMM in the pseudoscalar meson 
channel 





Fig. 1. Light-by-light contribution from intermediate pseu- 
doscalar meson exchanges. 



The light-by-light contribution due to exchanges of the 
light hadrons in the intermediate pseudoscalar channel to 
the muon AMM is shown in Fig. [T] and can be written in 
the form 1261 



00 00 



,LbL,PS 



x E 

a— 7T° ,77,77' 



J dq 2 J dq 2 J dtVl~t 2 \x 

00-1 
, F a . 7 . 7 . (g|;q 2 ,g 2 )F a » 7 . 7 (gf;g|,0) 
ql + M 2 a 



F a * 7 * 7 * {qi', <?i, 9I) F a * 7 * 7 



h 



el 



(20) 



where q-j = — (qi + q 2 ) ■ The functions I\ and I 2 are 
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h 



+ 



1 

D~i 



(9192) 



? ?(i-< s ))+(gf 



1 



Dx-D 7 



h = 2q 2 2 



(9192) 



(21) 



1 



(«? + (M 2 ) + M;(i-t 2 )) 

91(9192) - (9192), 



(22) 



are obtained [26] after averaging over the directions of 
muon momentum p 



1 

2^2 



(23) 



In (HU) and (HH) the notations arc (D 1 = (p + qif + M%, 

D2 = (p- 92) 



Ml) 



1 



Ri - 1 
IP"' 



1 



R2-I 
2MI 



1 



£>i-£> 2 / MaiftHgala; 

P9i\ , JI--R2) 2 
W = - {qiq2) ^Mj- 



arctan 



1 - zt 



(24) 



P92 

Di 

t = 



= (9192) 
(9192) 



8M2 



I9i| |92 
9i 92 



as = v 7 ! - < 2 , i2i = Wl 



4M,2 



4M, 2 



(1 - ife) (1 - R 2 ) , 



and is the muon mass (p 2 — —Mz) ■ 

It is instructive to investigate the pion contribution 
to the muon AMM o£ bL ' ,r for the 577(2) reduction of 
the nonlocal model (JXJ) - In this case, in the isospin limit, 
there are three model parameters: the current m c .„ and 
dynamical m^„ quark masses and the nonlocality param- 
eter A. Varying one parameter, say m^, in some region, 
one can fix other parameters by using as input the pion 
mass and the two-photon decay constant of the neutral 
piorfl. The pion mass is known with high accuracy, but 



4 The properties of the charged pion, namely the mass and 
the weak decay constant, are used as a rule for fixing the model 
parameters. In general, the difference between fixing the pa- 
rameter space for the charged pion or for the neutral pion is 




200 225 



250 



275 300 

m d , u [MeV] 



325 350 



Fig. 2. Light-by-light contribution to the muon AMM from 
the neutral pion exchange as a function of the dynamical quark 
mass. The band between dotted lines corresponds to the error 
interval for the pion two-photon width. 



the two-photon decay constant has an experimental un- 
certainty at the level of 3%. We investigate the depen- 
dence of on this uncertainty. Namely, the dynam- 
ical quark mass is taken in the typical interval 200-350 
MeV and then other parameters are fitted by the pion 
mass and the two-photon decay constant in correspon- 
dence with the pion lifetime given within the error range 
in [57]. The dependency of a^ bL,7r as a function of m^ u 
within the error interval for the two-photon decay con- 
stant is shown in Fig. [5] and, thus, we get our conservative 
estimate o£ bL ' ir ° = (5.01 ± 0.37) • 10~ 10 . 

The results for the nonlocal SU(3) model (JXJ) with the 
parameter sets G/_/y are given in Table [3] Note, that the 
results for a^ hL ' v are in the error range consistent with 
the results obtained in the nonlocal SU(2) model and that 
fl LbL,r; an( j qLIjL.jj arc onc or j cr i css than the contribution 

to the muon AMM due to the pion. 



set 


7T U 


V 


v' 


77 + 77' 


7T U + 77 + 77' 


Gi 


5.05 


0.55 


0.27 


0.82 


5.87 


Gii 


5.05 


0.59 


0.48 


1.08 


6.13 


Gin 


5.05 


0.53 


0.18 


0.71 


5.76 


Giv 


5.10 


0.49 


0.25 


0.74 


5.84 



Table 3. The contribution of pseudoscalar mesons to the muon 
AMM a bbL obtained within NAQM for different parametriza- 
tions. All numbers are given in 10~ 10 . 



due to small isotopic invariance breaking corrections induced 
by the electro-magnetic interaction and inequality of the non- 
strange quark masses. We believe that for the problem under 
consideration it is better to fix the model parameters by the 
neutral pion properties. 
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0.30 




2.0 3.0 

Q 2 [GeV 2 ] 



Fig. 3. The 7*7 — > ir° transition form factor in NxQM (this 
work, parameter set Gi is used), LMD+V [7J and NJL [12] 
models in comparison with the CELLO [28], CLEO [29] and 
BABAR [30] data. 



5 Constraints on F77 vertex and LbL 
amplitude 



A first phenomenological constraint on the anomalous ver- 
tex Fp» 7 » 7 * fgf; qf, gf) is put when both photons and the 
meson are on-shell. The vertex Fp 77 (— Mp;0, 0) is nor- 
malized by the experimentally measured two-photon de- 
cay widths 



p 77 



-Mp; 0,0) 



: 5p 77 



/ 647rr (P -> 77) 
(47ra) 2 M% 



A second phenomenological constraint is, that the ver- 
tex for special kinematics, when the meson and one of 
the photon are on-shell, Fp 77 - (— Mp;0, q 2 ), has to fit 
the data on the pseudoscalar meson transition form fac- 
tors available from the measurements of the CELLO [28] . 
CLEO [25J and BABAR [3DJISI] collaboration^. The tran- 
sition form-factors for tt, ry, r\' mesons evaluated in the 
NxQM for the G/ parameter set are shown in Figs. 13 1415 1 
One can see the reasonable agreement of the NxQM (this 
work) and LMD+V (VMDJJ [7] models predictions with 
existing experimental data, while the NJL model with the 
parameter set taken in [12] is not perfect in explanation 
of the data. 

A first QCD constraint [35] is, that at large photon 
virtualities q 2 + q 2 — > 00 the meson transition form factor 
has the asymptotic behavior as 



p 7 * 7 » 



-M 



As was discussed in [321133 1134] . new, very interesting data of 
the BABAR collaboration 30,31 at high momentum transfer 
of the photon do not affect very much the low energy observ- 
ables. 

6 Some details about the LMD+V and VMD models are 
given in Sect. 6. 



+ 
o 

CM 

0? 

c-i p. 



CELLO 
CLEO - 
BABAR 
NxQM - 
VMD 




2.0 3.0 

Q 2 [GeV 2 ] 

Fig. 4. The 7*7 — > r\ transition form factor in NxQM (this 
work, Gi set) and VMD [7J models in comparison with the 
CELLO [28] . CLEO [29] and BABAR [31] data. 



CELLO 
BABAR 
N X QM - 
VMD 





0.40 r 


I- 


0.35 \ 








0.30 - 








0.25 








0.20 








0.15 - 








0.10 


^_ 


0.05 - 




0.00 



0.0 




2.0 3.0 

Q 2 [GeV 2 ] 



Fig. 5. The 7*7 — > r\ transition form factor in NxQM (this 
work, Gi set) and VMD [7J models in comparison with the 
CELLO [21] and BABAR [3TJ data points. 



and for symmetric kinematics it fixes the asymptotic co- 
efficient. In the case of pion one has 



q 2 3 



6 7T 



+ 0(a 2 s )+0(q- 2 ) 



(25) 



Note, that even, if both in the NxQM and the LMD+V 
approaches the 1/q 2 behavior is satisfied and power cor- 
rections are taken into account, these models are still in- 
complete to describe the radiative corrections O (a s ) cal- 
culated within the perturbative QCD [315]. 

A second QCD constraint is to match the total light- 
by-light amplitude with three virtual {qui — 1,2,3) and 
one external magnetic (54 — > 0) field 

M {qlqlql) = a 2 N c Tr [q 4 ] A {q\, q\, q 2 ) 

to the perturbative QCD asymptotics for the configura- 
tions when some of photon virtualities are large [5] . There 
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are only two distinct kinematic regimes in the light-by- 
light scattering amplitudes: the Euclidean momenta of the 
three photons are comparable in magnitude q\ rs q 2 rs 
q\ 3> vIqcd i or one °f the momenta is much smaller than 
the other two q\ rj g| ^> g| 3> v4q Cd . The second limit 
was analyzed in [5] using the OPE with the result that the 
amplitude is factorized in the short-distance factor ~ 

times the nondiagonal correlator f VvlVj between the ax- 
ial and electromagnetic current in the external magnetic 
field 



Aps {(ft, 



(26) 



E 

a=3,8,0 



two arguments and in specific kinematics it is reduced to 



G 2 (8,8,8) = —j, 



3s^ 



G 2 (83,8,8) 



(2s- 



s 3 (33 - 4s) 
4s - s 3 + 2s (s - s 3 ) J (33,8, s) 
4 



S3 )ln^ + 
s 



G 2 (s, s, si) 



(4s 2 - 



In- 



si 



ssi (si - 4s) 
2s (4s - si) - 2s (s 2 + 2s 2 - 3ss x ) J (s, s, si) 



(29) 
(30) 

(31) 



where 



G 



(a) 



93>92.9l, 

1 



2 1 



W {3) = - W {s) = 
4' 



1 

12' 



W (0) = 



3' 



,(<») 



(gf ) is the longitudinal part of the VAV correlator 
with the corresponding flavor structure of the axial cur- 
rent, jf v — q^e" — q^e^ are the field strength tensors, 
the braces denote traces of the products of the matrices 
f? v . In (|2"tJ|) . the amplitude Aps is the part of the to- 
tal amplitude that is relevant to hadronic exchange in the 
pseudoscalar channel. The properties of the VAV correla- 
tor was studied in detail in [37J where it was shown that 
in the chiral limit of massless current quarks, due to prop- 
erties of the axial anomaly, in the perturbative QCD one 
has for all a 



(a) 



(27) 



moreover, these expressions for wf' 8 ^ (q 2 ) are exact QCD 
results for nonsinglet axial currents valid at any q 2 . This 
result was also confirmed within the effective instanton 
liquid model (N%QM) [35]. In [33] it was also demon- 
strated how the anomalous longitudinal part of the cor- 



relator 



,(°) 



q 2 ) is modified at low momenta (q 2 < ^ 2 ') 
due to the presence of the V a (1) anomaly. Thus within 
OPE one has, independently of a, the QCD asymptotic 
constraint 



where s$ = q 2 . The triangle function J (S3, s 2 , si) (see, for 
example, [40] ) is symmetric in all its arguments and for 
the kinematics considered it is equal to 



J(s3,s,s)= (32) 

21n 2 (A^)-4Li 2 (A T A)-ln 2 ( y ) + ^ , 
if 4s < s 3 , 



j_ 1 

S3 A 





n ■ t 


( dthx 


2 sin — 


la 


2 



where x = f (0 < x < 4) , y = ^ (0 < y < 1/4) , A = 
y/1 — 4y, Cl 2 (8) is the Clausen integral function 



Cl 2 (8) 



and Li 2 (77) is the dilogarithm (Spence) function 



Li a (»?) = - / dt \ 
Jo 1 



Asymptotically, when some photon virtualities are large, 
one has 



ri I 2 2 2 
G 2 (q 3 ,q 2 ,qi 



2A 9?«92»93>^QCD 



2 2 ' 



(28) 



The function G 2 (gf, gf , g 2 ) was given in [5J as a result 
of explicit calculations in perturbative theory for the light- 
by-light scattering amplitude where the photon-photon in- 
teraction is mediated by the loop of massless quarks. This 
function is symmetric under the permutation of the last 



G 2 (s 3 ,s,s) = 

ss 3 



o 



G 2 (s,s,sx) 



(33) 



(34) 
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and 



of the pion-exchange part of the total LbL amplitude as 



(*2{S3,S,S) = j 

S 3 



In 



«3 



— 2 In — + 10 In — +4+ -7T 




«3 



/ \ s<si 4 

G 2 {s,s,si) = 



•S3 



In [ — 

si 



(35) 




(36) 



These asymptotics have to play important role in con- 
straining the effective nonperturbative models when cal- 
culating the loops of dynamical quarks |13j . 

In some recent works, however, it was claimed and used 
a so-called "new QCD constraint" on the pion exchange 
LbL contribution to the muon g — 2 [^I fT^lfTD] . This " con- 
straint" is an attempt to restrict the behavior of the pion- 
two-photon vertex when the pion is off-shell. In particular, 
in |H| it was claimed that it would be followed from QCD 
that 



? oLMDV 



<zV,o: 



1 »^c 



(37) 



where /o is the pion-decay constant in the chiral limit 
and x is the quark condensate magnetic susceptibility in 
the presence of a constant external electromagnetic field. 
Let us show that (|3T|) is inconsistent with OPE and QCD 
constraints (|2"5 ]) -([M ]) , 

Indeed, in [5] the off-shell pion form factor is con- 
structed as (off-shell LMD + V model) 



poLMDV ( 2. 2 2\ 
_ fir -<£<& (ll 



(38) 



4 



4) 



3 (q 2 + M 2 ) (q\ + M 2 ) (q 2 + M?) (q 2 + Mf ) ' 

Ph {QhQiQl) = hi (ql + ql) 2 + h 2 q\ql+ 
+ h 3 (ql + ql) ql + h±q\ - h 5 (ql + ql) - h e ql + h 7 , 



where M(i 2 ) are the the masses of the lowest vector me- 
son states, and the coefficients in (l38l) are fixed from the 
above mentioned constraints. Specifically, h-j is due to 
the pion two-photon decay width, h\ = by require- 
ment ([25]), h 5 is from the fit of the CLEO data [29] . 
h 2 is related to higher power corrections in (|25[) and fi- 
nally (ft. 3 + h^) M^~ 2 M 2 ~ 2 = x is fixed from the "new con- 
straint" (|3"T)) . As a result from (|38p one get the asymptotics 



ot^ v (lUUl) 



goLMDV 



J 2,PS 



oLMDV 
2, PS 



(s,s,s) 
(s 3 ,s,s) 



q§ + A 

»»4cD 1 

~3 



/fo- 



ssa 



(39) 

(40) 
(41) 



Comparing the asymptotic coefficients in (|40p and (14"TT) 
with corresponding QCD coefficients in (|2"9l and ([33]) it is 
clear that the results based on the "constraint" ((37)) are 
in contradiction with the QCD asymptotics obtained in a 
model independent way. We have already commented in 
[13] on the origin of this confusion as an incorrect identifi- 
cation of the (PVV) correlator with (irVV) . In the latter 
case, it appears the physical hadronic pseudoscalar cur- 
rent that leads to much higher suppression (|42l) than it 
is shown in (|37[) where the local pseudoscalar current is 
used. The behavior (|37p is in variance with the prediction 
of the NxQM where from (ITS)) one has at large q 2 



F™>W,0) " 



exp 1 



f/A 2 ) 



(42) 



that does not violate the OPE results. 

At this point, for completeness, we remind that the 
LMD+V model [7] corresponds to the expression (13"81 at 
q\ = 0. The VMD model, used in particular in [I] to de- 
scribe the form factors of the 77 and rj' mesons, is defined 
as 



pVMD 



(o-,qi 



il) 



5P77^M 



A 2 m) 



(Q 2 2 



A 2 mY 



(43) 



where the parameters Am are taken from the fit of CLEO 



data [221: A„ 



776 ± 22 MeV, A„ = 774 ± 29 MeV, 



859 ± 28 MeV. 



As a final remark we note, that from the results of [371 
8 , 38 it is easy to get the correct expression that includes 
the magnetic susceptibility x in the asymptotic expansion 
of the total light-by-light amplitude. Indeed, by using the 
OPE as it was shown in [37] (see also [38]) due to sponta- 
neous breaking of the chiral symmetry there is the power 
correction to wl (q 2 ) that is linear in the current quark 
mass 



Aw L (q 2 ) 



16tt 2 m c \(qq)\x 



(44) 



Then, we have the power correction to the total light-by- 
light amplitude 



AG 2 (q 2 ,ql,q% 



9?~92>93»^qcd 32tt 2 m c \(gg)\x 

3 q 2 4 



(45) 



which is suppressed by the current quark mass and by 
extra power of large momentum squared. 

Thus our conclusion is that the effective models, in 
particular the instanton liquid based NxQM model, are 
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able to satisfy to most of phenomenological and QCD 
constraints and take into account the full kinematic de- 
pendence on virtualities of the mesons and photons. As to 
VMD based models, it is the LMD+V model [7] for the 
pion form factor that satisfies to many constraints. The 
most serious problem of the VMD models is that they do 
not take into account the dependence of the meson-two- 
photon vertex on the meson virtuality and thus the mo- 
mentum conservation is violated in this vertex. We discuss 
this point with more details in the next section. 



6 Results and comparison with other models 

The results of different approaches in calculation of the 
pseudoscalar meson contributions to the muon AMM from 
light-by-light scattering mechanism are given in Table |H 
For the pion contribution to the muon AMM in NxQM we 
take the numbers from the SU(2) version of the model. For 
77 and 77' we use the average of the results obtained with 
different parameter sets, see Table [3] 

One can conclude that within the NxQM the pseu- 
doscalar meson contributions to muon AMM are system- 
atically lower then the results obtained in the other works. 
This effect can be understood considering more carefully 
the off-shell meson-two-photon vertex Fp* 7 * 7 (p 2 ; p 2 , 0) en- 
tering the light- by-light contribution Eq. (f!?0")) . The com- 
parison of the NxQM predictions for the pion in this kine- 
matics (fl8l) with other models is presented in Fig. [6j One 
can see that the NxQM leads to stronger suppression of 
the vertex at all momenta. This explains why the results 
for a^ bL ' 7r obtained within the NxQM is smaller compar- 
ing with other calculations. 

The situation with 77's is more dramatic. For illustra- 
tion in Fig. [7J we present on the same plot the vertex 
Fp» 77 (p 2 ; 0, 0) in the timelike region p 2 < and the ver- 
tex Fp* 1 * 1 (p 2 ;p 2 , 0) in the spacelike region p 2 > as they 
look in the NxQM and VMD models. These two special 
kinematics match at zero virtuality p 2 = 0. The remark- 
able feature of this construction is that the first kinematics 
is connected with the decay of pseudoscalar mesons into 
two photons at physical points Fp 77 (— M 2 ^; 0, 0) = gp 77 , 
while the second kinematics is relevant for the light-by- 
light contribution to the muon AMM ([2D]l. Thus, the part 
of Fig.[7Jat p 2 < describes the transition of the pion-two- 
photon vertex from the physical points of meson masses 
to the point with zero virtuality, which is the edge point 
of the interval where the integrand of (f^Uj) is defined. In 
VMD type of models, including LMD+V model, there is 
no such dependence on the meson virtuality. Thus, the 
value of this vertex at zero meson virtuality is the same 
as the value of the vertex at the physical points of meson 



masses. 



= -Mi 



0,0) = 2^(0; 0,0). How- 
ever, the 77 and 77' mesons are much heavier than the pion 
and such extrapolation is too crude. One can see that for 
77 and particularly for 77' the difference between the values 
of the vertex at physical and zero virtuality points is large 
in the NxQM, F Pll {p 2 = -M^jO.O) > i^p 77 (0; 0, 0). 
Thus, the contributions of the 77 and 77' mesons to the 



muon AMM evaluated in NxQM are strongly suppressed 
as compared with the VMD results that can only be con- 
sidered as upper estimates of these contribution^. 




2.0 3.0 

p 2 [GeV 2 ] 

Fig. 6. The photon-pion vertex for the special kinematics when 
the pion and one photon are off-shell and have the same mo- 
menta, and another photon is real. The NxQM, Eq. (JTHJ) , is 
given by the solid line, the VMD, Eq. (|43[1 . is by the dashed 
line, the NJL result [12] is by dash-dotted line, and oLMD+V, 
Eq. (f38)l . is by the short dashed line. 
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p 2 [GeV 2 ] 

Fig. 7. Plots of the 7r°, 77 and rj vertices -Fp» 77 (p 2 ; 0, 0) in the 
timelike region and _Fp» 7 * 7 (p 2 ;p 2 , 0) in the spacelike region in 
NxQM model (thick lines) and VMD, Eq. (g3J, model (thin 
lines). The points with error bars correspond to the physical 
points of the meson decays into two photons. The VMD curves 
for 7T° and 77 are almost indistinguishable. 



7 Recently another attempt to take into account the full kine- 
matic dependence of the pion-photon vertex is given in the DSE 
approach 15 . In this paper, the meson's contribution is in- 
vestigated using the pion-pole approximation with an off-shell 
prescription. Although this prescription leads to a small sup- 
pression a full T-matrix calculation therein is highly desirable. 
We thank Richard Williams for clarifying this point. 
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TV 


V 




7T +?7 + 77 


\7"|\/l"rt & 
VIVID [D| 


0.(4 


1 9 A 

1.04 


i in 
1.19 


o.z / ^U.047 


fn tt rm 








o.u ^1 .07 


LMD+V, VMD [7] 


5.8(1.0) 


1.3(0.1) 


1.2(0.1) 


8.3(1.2) 


NJL [H] 


8.18(1.65) 


0.56(0.13) 


0.80(0.17) 


9.55(1.66) 


(LMD+V)',VMDl8] 


7.97 


1.8 


1.8 


11.6(1.0) 


N X QM [S] 


6.5(0.2) 








HM [16] 


6.9 


2.7 


1.1 


10.7 


DIP, VMD [10] 


6.54(0.25) 








DSE [J5] 


5.75(0.69) 


1.36(0.30) 


0.96(0.21) 


8.07(1.20) 


This work (NxQM) 


5.01(0.37) 


0.54 


0.30 


5.85 



Table 4. The light-by-light contribution to the muon AMM from the hadronic exchanges in the neutral pseudoscalar channel 
aj; L obtained in different works. All numbers are given in 10 -10 . Note, that as it was checked in [411113] . the (LMD+L)' model 
[S] numerically predicts 7.97 instead of 7.65 quoted originally in [BJ. 



7 Conclusions 



In this paper we have analyzed the pseudoscalar meson ex- 
change contributions of the light- by-light process to a^ S ' LbL 
in the nonlocal chiral quark model. The basic new ele- 
ments of our work are the inclusion of the full kinematic 
dependence of the pseudoscalar meson-photon vertex on 
virtualities of photons and mesons and the effect of cur- 
rent quark masses. This study extends the previous work 
[13] where the full kinematic dependence was considered 
for the pion in the chiral limit. The dependence of the ver- 
tices on the pion virtuality diminishes the result by about 
20-30% as compared to the case where this dependence is 
neglected. We find that the most significant effect occurs 
for the contributions of the 77 and 77' mesons to the muon 
AMM. In this case the results are reduced by factor about 
3 in comparison with the results obtained in other effec- 
tive quark models where the kinematic dependence was 
neglected. Thus our main conclusion is that within the re- 
alistic nonlocal chiral quark model the total contribution 
of pseudoscalar exchanges to 

a PS,LbL ifj 

by approximately 
factor 1.5 less than the most of previous estimates. 

We reviewed the phenomenological and QCD constraints 
on the hadronic light-by-light scattering amplitude. It is 
shown that the effective models like N%QM are able to 
fit most of these constraints, while the models based on 
the vector meson dominance violate the momentum con- 
servation law and lead to overestimation of the meson ex- 
changes contribution to the muon AMM. We have also 
demonstrated that one of the constraints for the LbL am- 
plitude recently discussed in the literature [5] is not justi- 
fied. 

An important next step in the investigations of the 
light-by-light scattering within the nonlocal chiral quark 
model is to perform an extension to the so-called com- 
plete calculation (see [42ll6lfT2]) that includes the scalar 
and axial-vector meson exchanges, as well as takes into 
account the quark and meson box diagrams. Due to con- 
tact terms arising in the nonlocal model, a calculation of 
the box diagrams is technically rather involved and will 
be presented elsewhere. 
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